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Abstract 



We report magnetotransport measurements in ballistic graphene/hexagonal boron nitride meso- 
scopic wires where the charge carrier mean free path is comparable to wire width W. Magne- 
toresistance curves show characteristic peak structures where the peak field scales with the ratio 
of cyclotron radius R c and wire width W as W/R c = 0.9 ± 0.1, due to diffusive boundary scat- 
tering. The obtained proportionality constant between R c and W differs from that of a classical 
semiconductor 2D electron system where W/R c = 0.55. 
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13 The presence of unique relativistic charge carriers, massless Dirac fermions, in monolayer 

14 graphene [JJ has lead to the emergence of novel transport phenomena such as half-integer 

15 [21 E] and fractional quantum Hall effect [H [5]. When the elastic mean free path of Dirac 

16 fermions becomes comparable to the sample size, the charge carrier transport mechanism 

17 enters a quasi-ballistic regime [6], and Dirac fermions exhibit more intriguing relativistic 
is transport phenomena such as Klein tunneling [7], evanescent wave transport [8], and sub- 

19 Poissonian shot noise [9J. In quasi-ballistic graphene devices, the probability of charge 

20 carrier scattering at sample boundaries becomes significantly larger than in the bulk region, 

21 because the probabilities of charge carriers scattering from charged impurities [JO], defects 

22 [IT] , longitudinal acoustic phonons (T2J [13], substrate optical phonons [l2~| HI], and flexural 

23 phonons [J_5] are significantly suppressed. To understand the fundamental transport proper- 

24 ties of quasi-ballistic graphene, it is highly desirable to characterize the effects of boundary 

25 scattering. 

26 In the mesoscopic wire system made from conventional semiconductor-based two-dimensional 

27 electron systems, it has been established that, if the charge carriers travel ballistically in 

28 the bulk region and are scattered diffusively at the sample boundary, the device exhibits 

29 anomalous magnetoresistance peaks owing to the magnetic commensurability effect between 

30 cyclotron radius R c and wire width W [HI US] • Moreover, since the anomalous magnetore- 

31 sistance peaks are due to diffuse boundary scattering, the probability of specular scattering, 

32 namely specularity parameter p, can be extracted by measuring anomalous magnetoresis- 

33 tance peaks [T61H9] . 

34 Recent development of a transfer technique of graphene on hexagonal boron nitride (h- 

35 BN) J20J [21] have made a breakthrough in the studies of ballistic transport phenomena 

36 in graphene [22]- Since the surface of h-BN is atomically flat and h-BN has high optical 

37 phonon energy as compared to SiC>2, graphene on h-BN exhibits extraordinary high mobility 

38 /i ~ 100,000 cm 2 /Vs and a long mean free path / m f p = (h/2e)/j^/n/7t ~ 1 /zm [221 123]. In 

39 conventional graphene on SiC>2 devices, fi is relatively low, and Z m f p is limited to 100 nm [T2] . 

40 Suspended graphene devices show high /i = 200, 000 cm 2 /Vs [2H I2S], but the observation of 

41 boundary scattering effect has not been reported [2S] • 

42 In this letter, we report on magnetotransport measurements in ballistic graphene/hexagonal 

43 boron nitride mesoscopic wires, in which the carrier mean free path is comparable to the 

44 wire width. We observed anomalous magnetoresistance curves with characteristic peak 
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45 structures where the peak field scales with the ratio of cyclotron radius R c and wire width 

46 W as W/R c = 0.9 ±0.1, which indicates the detection of diffusive charge carrier scattering at 

47 the graphene boundary. The obtained proportionality constant between R c and W contrasts 

48 that of the classical semiconductor two-dimensional electron system, in which W/R c = 0.55. 

49 In addition, from the analysis of magnetoresistance amplitude, nearly zero specularity at 

50 the graphene boundary is suggested. 

51 The graphene /hexagonal boron nitride (GBN) mesoscopic wire system was fabricated 

52 using the following lithography steps [201 I2TH29"] . First, relatively thick h-BN crystals (~ 10 

53 nm) were deposited on a Si wafer using a mechanical exfoliation technique. A monolayer 

54 graphene flake was deposited on a spin-coated polymethylmethacrylate (PMMA) layer and 

55 transferred on a h-BN crystal using an alignment technique under an optical microscope [20] . 

56 Monolayer thickness of a graphene flake was verified by Raman spectroscopy. Bar-shaped 

57 geometry was defined using standard electron-beam lithography and subsequent oxygen 

58 plasma etching. The electrical contacts were defined by electron-beam lithography followed 

59 by the evaporation of Pd (80 nm) and lift-off. The resist residues were removed by annealing 

60 at 300 °C in Ar/H 2 (97:3) gas flow for 6 h. 

61 In order to study the dependence of the transport properties on channel length L, we 

62 fabricated three devices for L = 0.6, 1.4, and 2.3 fim with channel width W = 1.0 /zm in a 

63 single graphene flake [inset in Fig. 1(a)]. Transport measurements were carried out using 

64 the standard lock-in technique with a small alternating current of 7 ac = 100 nA in a variable 

65 temperature insert at T = 1.5 ~ 300 K. A heavily doped Si substrate was used as a global 
ee back gate V g to tune carrier density n = C g (V g — Vbi rac ), where C g = 1.07 x 10~ 4 F/m 2 is the 
67 gate capacitance and VDi rac is the value of V g at the charge neutrality point [2j . A magnetic 
es field was applied perpendicularly to the sample surface. 

69 Figure 1(a) shows conductivity a as a function of for the L = 2.3 /xm device measured 

70 at T = 4, 16, 60, 120, 180, 240, and 300 K (from top to bottom). To extract conductivity a 

71 from the measured resistance, the contribution of contact resistance were eliminated using 

72 transfer length method [supplementary information]. Typical V-shaped dependence of a 

73 on V g was observed for all temperatures. Minimum conductivity was located close to zero- 

74 gate bias voltages Voirac ~ —0.3 V, indicating small unintentional doping (~ 10 10 cm -2 ). 

75 The carrier mobility limited by long-range scattering (described below) on hole side reached 

76 /i\ = 63, 000 cm 2 /Vs at T = 4 K. In the following discussions, we focus on the hole side of 
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77 the a because the hole mobility was higher than the electron mobility. 

78 To extract the mean free path, we fitted the conductivity using the transport formula 

79 and assuming long- and short-range scattering as the major contributors to charge carrier 
bo scattering [30, 3 lj , a -1 = (nep\)~ l + p s , where p\ is the mobility from long-range scattering, 
si and p s is the resistivity from short-range scattering, which are both independent of charge 

82 carrier density n [20]. In —35 V < V g < V, the experimental data are well fitted to the 

83 calculated conductivity [dashed curve in Fig. 1(a)]. In —50 V < V g < —35 V, the calculated 

84 conductivity deviates from the experimental data. This deviation can indicate the presence 

85 of an additional scattering mechanism, which is neither long- nor short-range scattering. 

ae From the obtained values of p\ and p s for varying temperatures T, we evaluated the 
87 mean free path for long-range scattering Zi ong = fi^fwKp\e~ x and the mean free path for 
ss short-range scattering / short = h^fnrr^ p^ 1 e~ 2 [21]. The values of /i ong and / s h or t as a function 

89 of V g for T = 4, 16,60, 120, 180,240, and 300 K are shown in Figs. 1(b) and (c) (from top 

90 to bottom). The value of Zi ong increased with |Vg| [Fig. 1(b)], whereas the value of l s hort 

91 decreased with \ V g \ [Fig. 1(c)]. As the temperature increased, both /i ong and / s hort decreased. 

92 These observations were qualitatively consistent with the previous transport measurements 

93 in graphene on Si0 2 . Quantitatively, however, the obtained values of l\ ong ~ 1.2 pm and 
a* 4hort ~ 2 /iin for Vg = —45 V were substantially larger than the typical values in graphene on 

95 Si02, which were around hundred-nanometer range [12] ; moreover, these values were larger 

96 than sample width W. 

97 Figure 2(a) shows two-terminal resistance R as a function of magnetic field B at V g = —45 

98 V and T = 4 K for the L = 2.3 pm and W = 1.0 pm device. As we increased B from zero, R 

99 increased with B and exhibited a maximum at -B max = ±0.21 T. As B was further increased, 

100 R decreased with B and exhibited a minimum at B m [ n = ±0.51 T. For larger B, R oscillated 

101 as a function of B, which can be attributed to Shubnikov-de Haas oscillation. 

102 In standard theory for ballistic transport phenomena in two-dimensional electron systems 

103 in a small magnetic field [B], if the cyclotron radius is larger than the sample width, R c > W, 

104 and if the electrons are scattered diffusively at the sample boundary, the boundary scattering 

105 effect leads to an increase in the backscattering probability as schematically shown by solid 
loe curves in the inset of Fig. 2(a). When the magnetic field is increased, the cyclotron diameter 
107 becomes smaller than the sample width 2R C < W. Under this condition, a skipping orbit of 
los charge carriers at the sample boundary is formed, and the backscattering probability is sup- 
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109 pressed as schematically shown by the dashed curves in the inset of Fig. 2(a). Moreover, for 
no this condition, the quantization of the cyclotron orbit, i.e., Landau quantization, is allowed 
in because 2R C < W, thus Shubnikov-de Haas oscillation can be observed. The observation 
112 of anomalous peak structures and Shubnikov-de Haas oscillation in Fig. 2 is thus qualita- 
n3 tively consistent with the expectations from the ballistic transport with diffusive boundary 
H4 scattering in the mesoscopic wire system [32] . 

H5 We studied R as a function of B for shorter channel lengths of L = 1.4 and 0.6 /zm 
lie [Fig. 2(b)]. The amplitude of anomalous magnetoresistance curve, AR = -R max — -Rmin 
n7 [indicated by solid and dotted arrows in Fig. 2(a), respectively], decreased with decreasing 
us channel length L [Fig. 2(c)]. In standard theory for diffusive boundary scattering, the 
n9 amplitude of anomalous magnetoresistance peaks is expected to decrease with L, because the 

120 probability of back scattering decreases. Therefore, the observed decrease in the amplitude 

121 of magnetoresistance peaks for L = 1.4 and 0.6 /iin devices is also consistent with the 

122 expectations from boundary scattering in the mesoscopic wire system. From the observations 

123 in Fig. 2, it can be concluded that we have observed magnetoresistance peak structures due 

124 to diffusive boundary scattering in a ballistic graphene mesoscopic wire system. 

125 For quantitative understanding of the observed magnetoresistance curves, we discuss 

126 here the positions of magnetoresistance peaks B max . For diffusive boundary scattering, 

127 the value of B max scales with the ratio of cyclotron radius R c (B,n) to the channel width 

128 W as W/R c (B,n) = a, where a is the proportionality constant. For Dirac fermions in 

129 monolayer graphene, the cyclotron radius can be written as R c (B,n) = ^y^im, where e 



130 is an elementary charge and h is Planck's constant. Therefore, the expected positions of 

131 resistance peaks in graphene can be expressed as 



132 To study the relationship between R c and W, we measured magnetoresistance curves at 

133 varying gate-bias voltages for the sample with (L, W) = (2.3 /mi, 1.0 /im). Fig. 3(a) shows 

134 two-terminal resistance R as a function of magnetic field B at V g = —50, —48, • • • , —28 V 

135 (bottom to top). In this plot, each curve was offset vertically so that peak positions were 

136 separated by 0.01 kf2. At V g = —50 V (bottom curve), we observed the anomalous peak 

137 structures at -B max = ±0.21 T. As we increased gate-bias V g , the positions of peak structures 

138 were shifted to smaller B. At V g = —28 V (top curve), the peak positions were located at 




B, 




(1) 



= a — 



e W 
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139 B mSLX = ±0.14 T. The peak positions were compared with the calculated B max using equation 

140 (1). As shown by the colored areas in Fig. 3(a), the peak positions were well fitted with the 
mi calculated positions for a = 0.9 ± 0.1. 

142 To investigate the dependence of a on the channel width, we studied another sample for 

143 larger channel width (L,W) = (3.1 fim, 1.5 /mi) [Fig. 3(b)]. Compared with the device 

144 with smaller W [Fig. 3(a)], the values of -B max were decreased [Fig. 3(b)]. Moreover, the 

145 peak positions were also fitted by the equation (1) for a = 0.9 ± 0.1, as indicated by the 
we colored area in Fig. 3(b). These analysis were also valid for other samples for (L, W) = (1.4 
147 /mi, 1.0 fim) and (2.5 /xm, 0.95 /im). These observations indicate the independence of the 
us proportionality constant a on channel width W. 

149 The obtained proportionality constant a was larger than that of classical massive charge 

150 carriers in the conventional semiconductor two-dimensional electron system (2DES) where 

151 a = 0.55 [12]. For comparison, we plot the expected positions of peak structures for a = 0.55 

152 (blue curve) in Figs. 3(a) and 3(b). The curves for a = 0.55 could not explain the position 

153 of peak structures. The conventional proportionality constant between R c and W has to be 

154 modified to explain the observed transport phenomena in graphene. 

155 The differing value of Wj R c for graphene compared to that of conventional semiconductor 

156 2DES indicate that the scattering process of massless Dirac fermions at sample boundaries 

157 can be different from that of conventional 2DEGs. As a possible explanation for the differing 

158 values of W/R c , the interference of charge carriers at sample boundaries can be pointed out 

159 [33]. Since the charge carriers in graphene are chiral, when the charge carrier scattering 
wo at sample boundaries changes the momentum, the phase of their wave function are shifted, 
lei This phase shift of wave function at sample boundaries induces complex interference of 

162 charge carriers, which can alter the transport properties of ballistic graphene [33]. Moreover, 

163 the first-principle numerical calculations of the magnetoresistance in graphene nanoribbon 

164 suggested W/R c = 0.8 — 1.0 [31], which seems be consistent with our results. Albeit, for the 

165 exact understandings of the reason for the differing values of W/R c , more elaborated works 
we will be needed in the future. 

167 Finally, we study the magnetoresistance curves for varying temperatures T and gate-bias 

168 voltages V g . Fig. 4(a) shows the amplitude of magnetoresistance AR, long-range mean free 

169 path h ong , and short-range mean free path / s h or t as a function of T at V g — —35 V. As T was 

170 increased from 4 K to 240 K, the value of AR was decreased, but was retained to T = 240 
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171 K [Fig. 4(b)]. This observation implies that the channel region of graphene mesoscopic wire 

172 was still ballistic at T = 240 K, because of the relatively small electron-phonon scattering 

173 probability in monolayer graphene compared to that of conventional semiconductors [T2] , 

174 which demonstrate the characteristic electronic properties of graphene. 

175 Fig. 4(b) shows AR, l\ ong , and / s h or t as a function of at T — 4 K. The value of AR 

176 increased with |V^| from 20 to 40 V and saturated for 40 to 50 V [Fig. 4(b)]. In principle, 

177 AR increases with the mean free path in the bulk region, and AR saturates as the mean 

178 free path in the bulk region exceeds the sample width [6J. In Fig. 4(b), as we increased 

179 \V g \, /^ng increased, and for 40 V < \V g \ < 50 V, Zi ong exceeded the sample width [indicated 
ibo by the dotted line in Fig. 4(b)]. On the other hand, Z sn0 rt was decreased with \V g \ and 
ibi were larger than W. The standard theory for ballistic transport suggests that, if specularity 

182 parameter p is zero, AR emerges when the carrier mean free path in the bulk region becomes 

183 comparable to the sample width W |16| . If we assume that the charge carrier scattering in 

184 the bulk region is dominated by long-range scattering, the emergence of AR at /i on g ~ 1 /-tni 

185 indicates the almost zero p in our device. 

186 In summary, we conducted magnetotransport measurements in quasi-ballistic graphene/hexagonal 

187 boron nitride mesoscopic wires, in which the carrier mean free path is comparable to the wire 

188 width. Magnetoresistance curves show characteristic peak structures where the peak field 

189 scales with the ratio of cyclotron radius R c and electronic wire width W, which indicates the 

190 detection of a diffusive scattering effect at the sample boundary. The analysis suggests that 

191 the values of R c and W scale as R c /W = 0.9 ±0.1. This is in contrast to the semiconductor 

192 two-dimensional electron system where R c /W = 0.55. The analysis also suggests nearly zero 

193 probability for specular scattering at the sample boundary. These findings are fundamental 

194 step forward in understanding the effects of boundary scattering on the transport properties 

195 of nanostructured graphene devices such as graphene nanoribbons (35] and graphene single 
we electron transistors [36j |37] . 
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v g (V) v„ (V) 

FIG. 1. (Color online) (a) Conductivity a as a function of back-gate bias voltage V g measured 
at T = 4, 16, 60, 120, 180, 240, and 300 K (top to bottom). For clarity, each curve was offset by 2 
mS. Dashed curves are the calculated conductivity based on the Boltzmann equation discussed in 
the main text. The inset shows the atomic force microscopy image of the GBN sample studied in 
this work. The white region indicates the region of Pd metal electrodes. The white dashed lines 
indicate the outline of the graphene flake, (b) Short-range mean free path Z s hort and (c) long-range 
mean free path l\ ong as a function of gate-bias voltage V g , for T = 4, 16, 60, 120, 180, 240, and 300 
K (top to bottom). 
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FIG. 2. (a) Two-terminal resistance R as a function of B measured at T = 4 K for the device with 
channel length L = 2.3 (im with applied gate-bias voltage V s = —45.0 V. (inset) The schematic of 
electro-trajectories in mesoscopic wire for (solid curve) R c > W and (dotted curve) R c < W. (b) 
R as a function of B measured at T = 4 K for (top curve) L = 1.4 fim and (bottom curve) L = 0.6 
fim. (c) The amplitude of magnetoresistance signal Ai? vs. channel length L. 
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FIG. 3. (Color online) (a) -R as a function of B measured at T = (a) 4 and (b) 1.6 K for the 
device with (L, W) = (a) (2.3 fim, 1.0 ^m) and (b) (3.1 ^m, 1.5 /im) with applied gate-bias voltage 
V g = —50, —48, • • • , —28 V (bottom to top). Each curve was offset vertically so that peak positions 
were separated by 0.01 k^.The blue dotted curves and colored area indicate the expected peak 
positions for a = 0.55 and 0.9 ±0.1 respectively. 
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FIG. 4. (a) AR (top panel), Zi ong (blue curve in bottom panel), and Z s hort (black curve in bottom 
panel) as a function of T measured in the device with channel length L = 2.3 fim at V g = —35 V. 
(b) AR (top panel), l\ ong (blue curve in bottom panel), and Z s hort (black curve in bottom panel) as 
a function of V g at T = 4 K. The dotted line is the channel width of GBN mesoscopic wire. 
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